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Abstract 

The Skyrme model is a low-energy effective field theory for QCD, where the 
baryons emerge as soliton solutions. It is, however, not so easy within the stan- 
dard Skyrme model to reproduce the almost exact linear growth of the nuclear 
masses with the baryon number (topological charge), due to the lack of Bogo- 
molny solutions in this model, which has also hindered analytical progress. Here 
we identify a submodel within the Skyrme-type low energy effective action which 
does have a Bogomolny bound and exact Bogomolny solutions, and therefore, at 
least at the classical level, reproduces the nuclear masses by construction. Due 
to its high symmetry, this model qualitatively reproduces the main features of the 
liquid droplet model of nuclei. Finally, we discuss under which circumstances the 
proposed sextic term, which is of an essentially geometric and topological nature, 
can be expected to give a reasonable description of properties of nuclei. 

1 Introduction 

The Skyrme model G] is an effective low-energy action for QCD [2 |, where the pri- 
mary ingredients are meson fields, whereas baryons appear as solitonic excitations, and 
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the baryon number is identified with the topological charge. 
The original Skyrme Lagrangian has the following form 

L = L 2 + L 4 + L , (1) 

where 

f 2 

L 2 = - J -^ Tr (U^d„U U^d^U) (2) 

is the sigma model term, and a quartic term, referred as Skyrme term, has to be added 
to circumvent the standard Derrick argument for the non-existence of static solutions, 

U = Tr ([U*W, U^d^U} 2 ). (3) 

Here U is a 2 x 2 matrix- valued field with values in the group SU(2). The last term, 
which is optional from the point of view of the Derrick argument, is a potential 

L = -n 2 V(U,Ul), (4) 

which explicitly breaks the the chiral symmetry. Its particular form is usually adjusted 
to a concrete physical situation. The model has two constants, the pion decay constant 
/„• and the interaction parameter e. Additional constants may appear from the poten- 
tial. 

The modern point of view on the Skyrme model is to treat it as an expansion in deriva- 
tives of the true non-perturbative low-energy effective action of QCD, where higher 
terms in derivatives have been neglected. However, as extended (solitonic) solutions 
have regions where derivatives are not small, there is no reason for omitting such terms. 
Therefore, one should take into account also higher terms. In fact, many generalized 
Skyrme models have been investigated [3J, [4|, [5|, [6|, 

L = L 2 + L 4 + L + ..., (5) 

where dots denote higher derivatives terms. A simple and natural extension of the 
Skyrme model is the addition of sextic terms, among which one is rather special. 
Namely, we will consider the expression 

L 6 = ^ (Tr [er^rtdpU tfd v U U^d p U)f . (6) 

In standard phenomenology, the addition of this term to the effective action represents 
the inclusion of the interactions generated by the vector mesons u. In fact, this term 
effectively appears if one considers a massive vector field coupled to the chiral field via 
the baryon density [7|. Further, this term is at most quadratic in time derivatives (like 
the quartic Skyrme term) and allows for a standard time dynamics and hamiltonian 
formulation. In addition, it leads to a significant improvement in the Skyrme model 
phenomenology when applied to nucleons. Indeed, as explained first in Q, once the 
sextic term is present, it becomes the main responsible for stabilization, and then the 
quartic contribution changes sign as it corresponds to the scalar exchange it represents 
(solving an old puzzle). This compensation with the quadratic term holds also for the 
moments of inertia, when the rotation of all the mass is taken into account, as it should 
in the classical computation. 

In this letter we want to study the model restricted to the potential and the sextic term, 
Lqq = Lq + Lq, because this submodel has some unique properties. First of all, it has a 



2 



huge amount of symmetry [ 8 ] and, therefore, it is integrable in the sense of generalized 
integrability |9| (its symmetries and integrability properties shall be discussed in detail 
in a separate publication; its symmetries are also important for its rather close relation 
to the liquid droplet model of nuclei, as we shall discuss at length in the last section). As 
a consequence, the model has infinitely many exact solutions in all topological sectors, 
such that both energies and profiles can be determined exactly. Finally, the model has 
a Bogomolny bound which is saturated by all the exact solutions we construct below. 
The existence of static solutions which saturate a Bogomolny bound is very welcome 
for the description of nuclei, for the following reasons. Firstly, the resulting soliton 
energies obey an exactly linear relation with the baryon charges. Physical nuclei are 
well-known to obey this linear law with a rather high precision. Secondly, binding en- 
ergies of higher solitons are zero, again as a consequence of their Bogomolny nature. 
This conforms rather well with the binding energies of physical nuclei, which are usu- 
ally quite small (below the 1% level). Thirdly, the forces between sufficiently separated 
solitons are exactly zero. This result is a consequence of another crucial feature of our 
solitons, namely their compact nature. Again, this absence of interactions, although not 
exactly true, is a rather reasonable approximation for physical nuclei, given the very 
short range character of interactions between them. 

So we find a rather striking coincidence between some qualitative features of nuclei, 
on the one hand, and properties of our classical soliton solutions, on the other hand. 
One important question is, of course, whether this coincidence can be maintained at 
the quantum level. A detailed investigation of the quantization of the model is beyond 
the scope of this letter, but we shall comment further on it in the discussion section. In 
any case, the model seems to correspond to a rather non-trivial "lowest order" effective 
field theory approximation to nuclei which already reproduces some of their features 
quite well. We also want to remark that part of the pseudoscalar meson dynamics is 
possibly taken into account already by the potential Lo, which breaks the chiral sym- 
metry, as goldstone condensation. 

All the unique properties of the model may be ultimately traced back to the geometric 
properties of the proposed term Lq, i.e., to the fact that it is the square of the pullback 
of the volume form on the target space three-sphere S 3 (we remind that as a manifold 
SU(2) ~ S 3 ) or, equivalently, the square of the topological (baryon) current. We re- 
mark that models which are similar in some aspects, although with a different target 
space geometry, have been studied in IflOll . IfTTII . Further, the model studied in this 
letter, as well as its "baby Skyrme" version in 2+1 dimensions have already been intro- 
duced in 02) • There, the main aim was a study of more general properties of Skyrme 
models in any dimension. Concretely, the limiting behaviour of the full generalized 
Skyrme model for small couplings of the quadratic and quartic terms L2 and L4 was 
studied numerically. In addition, an exact solution for the simplest hedgehog ansatz 
was constructed, both in 2 and in 3 dimensions. For the three-dimensional solution, a 
rather complicated potential was chosen in order to have exponentially localized solu- 
tions, whereas in this letter we shall focus on the case of the simple standard Skyrme 
potential, which naturally leads to compact solitons. Besides, our main purpose is to 
make contact with the phenomenology of nuclei. The 2+1 dimensional baby Skyrme 
version of the model has been further investigated in iTPJll and recently in 1(1311 . with re- 
sults which are qualitatively similar to the ones we shall find in the sequel (e.g. compact 
solitons, infinitely many symmetries, Bogomolny bounds). 
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2 Exact solutions 



The lagrangian of the proposed restriction of the Skyrme model is 

L oe - ^ (Tr {tT^rfdJJ U^3 V U U^d p U)f - » 2 V(U, rf). (7) 

We start from the standard parametrization for U by a real scalar field £ and a three 
component unit vector n (r are the Pauli matrices), 

U = e**'*. 

The vector field may be related to a complex scalar u by the stereographic projection 

ft = - — \—r7T (u + u, -i(u - u), \u\ 2 - l) 
1 + \u\ z v ' 

giving finally (r± = (l/2)(n ± rr 2 )) 

where the SU(2) matrix field is 

W = (l + uu)~i 



1 iu 
iu 1 

and obviously cancels in the lagrangian. Using this parametrization we get (u p = d p u, 
etc.) 

A 2 sin 4 £ 

(i + H 2 ) 4 

where we also assumed that the potential only depends on tr U. The Euler-Lagrange 
equations read (V£ = c^U) 



L ™ = - n , L ,|" {^ Vp °t v u p u a f - i?V{0 (8) 



A2sin2? r ^(sin 2 ^^)- M 2 y e = o, 



(i + M 2 ) 4 

d - U + M 2 )V _ °' 

where 

„ _ d(e a ^^u p u a ) 2 _ d(e a ^^ u u p u a ) 2 

These objects obey the useful formulas 

H^u" = H p u» = 0, K p e = K p u» = 0, H p e = K p u" = 2(e av P^ v u p u a f. 

We are interested in static topologically non-trivial solutions. Thus u must cover the 
whole complex plane (n covers at least once S 2 ) and £ e [0, 7r] . The natural (hedgehog) 
ansatz is 

£ = £(r), u(6,4>)=g(6)e in +. 
Then, the field equation for u reads 

J^p. f 111 "\ 99j Q 



sin<9 \(l+g 2 ) 2 sm9J (1 + sin 2 9 
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and the solution with the right boundary condition is 

Q 

g{9) = tan-. 

Observe that this solution holds for all values of n. The equation for the real scalar 
field is 

n 2 X 2 sin 2 £ ( sin 2 jj r \ 2 

2r2 ^— j-A«V{ = 0. 

This equation can be simplified by introducing the new variable z = , 

sin 2 e^(sin 2 a,)-^=0, (9) 

and may be integrated to 

^sm^g = V, (10) 

where we chose vanishing integration constant to get finite energy solutions. Now, we 
have to specify a concrete potential. The most obvious choice is the standard Skyrme 
potential 

F=^Tr(l-L0 -> V(0 = l-cos£. (11) 

Thus, 

sin 2 £^ =±V2(l-cosO =>■ / ^-|- = ±2(z - z ). 



sin£/2 



o3 £ ,3, 



The general solution reads 
Imposing the boundary conditions for topologically non-trivial solutions we get 



cos 6 - =±-(z-z ). 



t = [ 2 arccos z e [0, f ] (U) 

\ z>|. 



The corresponding energy is 

1,2 ^4 



E = J d 3 x (~ (i + |^4 (Vr0 2 (VenV^ - V 0W V e S) 2 + M 2 ^ . (13) 



Inserting the solution for u and dTOb we find 

r /A 2 n 2 sin 4 ^ 
B - 47r/ r 2 dr — k 2 r+f x 2 V 



4r 4 

4tt • 2/i 2 / r 2 (W(£(r)) = 4\/27r^A|n| / dsV(f (*)) 



= 8%/2^A|n| jj 4/3 - (^y\ 2 Adz=^^fx\\n\. (14) 

The solution is of the compacton type, i.e., it has a finite support (compact solutions 
of a similar type in different versions of the baby Skyrme models have been found 
in |[T3l . lfT31 ). The function £ is continuous but its first derivative is not. The jump 
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of the derivative is, in fact, infinite at the compacton boundary z = 4/3, as the left 
derivative at this point tends to minus infinity. Nevertheless, the energy density and 
the topological charge density (baryon number density) are continuous functions at the 
compacton boundary, and the field equation (|9]i is well-defined there, so the solution is 
a strong solution. The reason is that £ 2 always appears in the combination sin 2 £ £ z , and 
this expression is finite (in fact, zero) at the compacton boundary. We could make the 
discontinuity disappear altogether by introducing a new variable £ instead of £ which 
satisfies £ 2 = sin 2 £ £ z . We prefer to work with £ just because this is the standard 
variable in the Skyrme model. 

In order to extract the energy density it is useful to rewrite the energy with the help of 
the rescaled radial coordinate 



y/2fx \ 3 / 3\n\z 



4A 



(15) 



like 

E = 8V2fiX ^ 47r ^' ' dff 2 (l - \n\-if 2 ) 

such that the energy density per unit volume (with the unit of length set by f) is 

£ = 8V2^A(1 - H~if 2 ). (16) 
In the same fashion we get for the topological charge (baryon number), see e.g. chapter 

i.4 of men 

1 f d 3 x - Sm *t ie mnl £ m u n iii (17) 



r2 



(i + H 2 ) 
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MS 



sign(n)— 4-7T / dff (1 — |n| 3 r ) 



o 

and for the topological charge density per unit volume 

B = sign(n)A;(l-\n\-ir 2 )z. (18) 

Both densities are, of course, zero outside the compacton radius f = |n| 3. We remark 
that the values of the densities at the center f = are independent of the topological 
charge B — n, whereas the radii grow like ni. For n = 1, we plot the two densities 
in Fig. CJ, where we normalize both densities (i.e., multiply them by a constant) such 
that their value at the center is one. 

We now want to compare the phenomenological parameters of our model (masses and 
radii) to the corresponding values for physical nuclei. One should keep in mind, of 
course, that the comparison is done at the purely classical level, and all quantum cor- 
rections are absent. First, observe that the energy of the solitons is proportional to the 
topological (baryon) charge 

E = E \B\, 
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Figure 1 : Normalized energy density (left figure) and topological charge density (right 
figure) as a function of the rescaled radius f, for topological charge n=l. For \n\ > 1, 
the height of the densities remains the same, whereas their radius grows like |n| s 



where Eq = 64v27T/xA/15. Such a linear dependence is a basic feature in nuclear 
physics. Let us fix the energy scale by setting Eq = 931.75 MeV. This is equivalent 
to the assumption that the mass of the B = 4 solution is equal to the mass of He 4 , 
which is usually assumed because the ground state of He 4 has zero spin and isospin. 
Therefore, corrections to the mass from spin-isospin interactions are absent. In table 
([D we compare the energies of the solitons in our model with the experimental values. 
We find that the maximal deviation in our model is only about 0.7%. For the numeri- 
cal determination of soliton masses in current versions of the Skyrme model we refer 
to ifTTl (the standard massive Skyrme model) and to lfl~8l (the vector Skyrme model, 
where a coupling of the Skyrme field to vector mesons is used instead of the quartic 
Skyrme term for the stabilisation of the Skyrmions). There, typically, the Skyrmions 
with low baryon number are heavier by a few percent, whereas they reproduce the linear 
growth of mass with baryon number for higher baryon number. In lfl9l the Skyrmion 
masses have been determined with the help of the rational map approximation |20| for 
Skyrmions, with similar results. 

Secondly, the sizes of the solitons can be easily computed and read 



Rb = -Royl-Bl, i?o = 



which again reproduces the well-known experimental relation. The numerical value is 
fixed by assuming Rq = 1.25 fm. 



3 Bogomolny bound 

Now, we show that our solitons are of the BPS type and saturate a Bogomolny bound. 
Let us mention here that a Bogomolny bound also exists for the original Skyrme model 
L2 + L4, but it is easy to prove that non-trivial solutions of this model cannot saturate 
the bound (see, e.g., [ 16]; this bound has been found already by Skyrme himself [ 1 1). 
The energy functional reads 

E= J ^(^^^^nutr+^v) = 
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-^06 


& experiment 


1 


931.75 


939 


2 


1863.5 


1876 


3 


2795.25 


2809 


4 


3727 


3727 


6 


5590.5 


5601 


8 


7454 


7455 


10 


9317.5 


9327 



Table 1: Energies of the soliton solutions in our model (Eqq), compared with the ex- 
perimental masses of physical nuclei. All numbers are in MeV 



,, 2/j.Xsm 2 £\/V _ 



± (2A/«r 2 



* t y d a; irT]# e ^ 

-i /" 3 sin 2 £^/V ? , 



2A^7r 2 Ci|B| (19) 



where B is the baryon number (topological charge) and the sign has to be chosen 
appropriately (upper sign for B > 0). If we replace W by one, then the result (i.e., 
the last equality in ( fT9b > follows immediately (and the constant C\ = 1). Indeed, for 
V = 1 the expression in brackets is just the topological charge ( fTTb . An equivalent 
derivation, which shall be useful below, starts with the observation that this expression 
is just the base space integral of the pullback of the volume form on the target space 
S 3 , normalized to one. Further, while the target space S 3 is covered once, the base 
space S 3 is covered B times, which implies the result. The same argument continues 
to hold with the factor \/V present (remember that V — V(£)), up to a constant C\. 
Indeed, we just have to introduce a new target space coordinate £ such that 



sm 



(VW)d^ = C lS m 2 ^dl (20) 



The constant C% and a second constant C%, which is provided by the integration of Eq. 
(l20l . are needed to impose the two conditions £(£ = 0) = and £(£ = tt) = it, which 
have to hold if £ is a good coordinate on the target space S 3 . Obviously, C\ depends 
on the potential V(£). Specifically, for the standard Skyrme potential V = 1 — cos£, 
C x is 

r , _ 32^2 

as may be checked easily by an elementary integration. We remark that an analogous 
Bogomolny bound in one lower dimension has been derived in J2T|, IfMl for the baby 
Skyrme model. 

The Bogomolny inequality is saturated by configurations obeying the first order Bogo- 
molny equation 



(1 + M 2 ) 2 

which, in the case of our ansatz, reduces to the square root of equation ( fTOb . The 
saturation of the energy-charge inequality by our solutions proves their stability. It is 
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not possible to find configurations with lesser energy in a sector with a fixed value of 
the baryon charge. 

4 Discussion 

In this letter we proposed an integrable limit of the full Skyrme model which consists 
of two terms: the square of the pullback of the target space volume (topological den- 
sity) and a non-derivative part, i.e., a potential. Both terms are needed to circumvent 
the Derrick argument. Then we explicitly solved the static model for a specific choice 
of the potential (the standard Skyrme potential). The resulting solitons satisfy, in fact, 
a Bogomolny equation. These exact Bogomolny solutions provide a linear relation 
between soliton energy (=nuclear mass) and topological charge (=baryon number B), 
reproducing the experimental nuclear masses with a high precision. Besides, these 
solitons have the remarkable property of being compact, which allows to define a strict 
value of the soliton size (=nuclear radius). The resulting radii R, too, follow the stan- 
dard experimental relation R <~ | | 1 / 3 with a high precision. 

These findings lead to the question of the nature and quality of the approximation which 
our model provides for the properties of physical nuclei. Obviously, the model as it 
stands cannot reproduce all features of nuclei, even qualitatively, because some essen- 
tial ingredients are still missing. First of all, the binding energy of higher nuclei is zero 
due to the Bogomolny nature of the solutions. Although not entirely correct for physi- 
cal nuclei, this is, however, not such a bad approximation, because the nuclear binding 
energies are known to be rather small. Their smallness is, in fact, one of the motiva- 
tions for the search for theories which saturate a Bogomolny bound. Secondly, there 
are no pionic excitations, because the corresponding term in the Lagrangian is absent. 
This is related to the complete absence of forces between separated, non-overlapping 
solitons. The absence of forces is a direct consequence of the compact nature of these 
solitons, because several non-overlapping solitons still represent an exact solution of 
the field equations. Physical nuclei are not strictly non-interacting, but given the very 
short range character of forces between nuclei, the absence of interactions in the model 
may, in fact, be welcome from a phenomenological point of view, within a certain 
approximation. Further, for physical nuclei a finite radius may be defined with good 
accuracy, so the compact nature of the solitons may be a virtue also from this point of 
view. 

For the energy density we find that it is of the core type (i.e., larger in the center and 
decreasing towards the boundary), see Fig. 1 . The baryon density profile is again of the 
core type, but flatter near the center, and with a smaller and more pronounced surface 
(=region where the density decreases significantly). For physical nuclei the density 
profile is quite flat (almost constant) and for some nuclei even with a shallow valley 
in the center, so here the phenomenological coincidence is reasonable but not perfect. 
Let us also mention that the independence of the profile heights of the baryon number 
conforms well with the known properties of nuclei. 

Our results for the profiles, however, have to be taken with some care. First of all, they 
depend on the form of the potential term, in contrast to the linear mass-charge relation 
(which holds for all potentials) or the compact nature of the solitons (which holds for 
a wide class of potentials). The second argument is related to the huge amount of sym- 
metry of the model. Indeed, for the energy functional for static field configurations, the 
volume-preserving diffeomorphisms on the three-dimensional base space are a subset 
of these symmetries. In physical terms, all deformations of solitons which correspond 
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to these volume-preserving diffeomorphisms may be performed without any cost in 
energy. But these deformations are exactly the allowed deformations for an ideal, in- 
compressible droplet of liquid where surface contributions to the energy are neglected. 
These symmetries are not symmetries of a physical nucleus. A physical nucleus has 
a definite shape, and deformations which change this shape cost energy. Neverthe- 
less, deformations which respect the local volume conservation (i.e., deformations of 
an ideal incompressible liquid) cost much less energy than volume-changing deforma- 
tions, as an immediate consequence of the liquid droplet model of nuclear matter. 
This last observation also further explains the nature of the approximation our model 
provides for physical nuclei. It reproduces some of the classical features of the nuclear 
liquid droplet model at least on a qualitative level, and the huge amount of symmetries 
of the model is crucial for this fact. Its soliton energies, e.g., correspond to the bulk 
(volume) contribution of the liquid droplet model, with the additional feature that the 
energies are quantized in terms of a topological charge. 

In other words, the model provides, besides a conceptual understanding with exact 
solutions, a new starting point or "zero order" approximation which is different from 
other approximations. It already covers some nuclear droplet properties of nuclear 
matter, and is topological in nature. For a more quantitative and phenomenological 
application to nuclei, obviously both the inclusion of additional terms and the quanti- 
zation of some degrees of freedom are necessary. 

So let us briefly discuss the question of possible generalizations of the model. A first, 
simple generalization consists in the choice of different potentials. The resulting soli- 
tons continue to saturate a Bogomolny bound, therefore the linear relation E ~ \B\ 
between energy and baryon number continues to hold. The energy and baryon charge 
densities for a spherically symmetric ansatz (hedgehog), and even the compact or non- 
compact nature of the solitons, however, will depend on the specific form of the poten- 
tial. 

A further generalization consists in including additional terms in the Lagrangian (like 
the terms L 2 and L4 of the standard Skyrme model) which we have neglected so far. 
From the effective field theory point of view there is no reason not to include them. If 
we omit, e.g., the kinetic term for the chiral fields, then there are no obvious pseudo- 
scalar degrees of freedom (77, n). These additional terms break the huge symmetry of 
the original model, such that the solitons now have fixed shapes. In order to describe 
nuclei, these shapes should be at least approximately spherically symmetric. A detailed 
investigation of this issue is beyond the scope of the present letter, but let us mention 
that at least under simple volume-preserving deformations from a spherical to an ellip- 
soidal shape both the E2 term and the E4 term energetically prefer the spherical shape. 
Further, the reasonable qualitative success of the restricted model might indicate that 
the additional terms should be small in some sense (e.g., their contribution to the total 
energy should not be too big). This opens the possibility of an approximate treatment, 
where the solitons of the restricted model L 6 provide the solutions to "zeroth order" 
(with all the topology and reasonable energies already present), whereas the additional 
terms provide corrections, which may be adapted to the experimental energies and 
shapes of nuclei. 

Further, a more realistic treatment certainly requires the investigation of the issue of 
quantization. We emphasize again that the rather good phenomenological properties of 
the model up to now are based exclusively on the classical solutions, and it is a different 
question whether quantum corrections are sufficiently small or well-behaved such that 
this success carries over to the quantized model. A first step in this direction consists in 
applying the rigid rotor quantization to the (iso-) rotational degrees of freedom, as has 
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been done already for the standard Skyrme theory [22], for some recent applications 
to the spectroscopy of nuclei see e.g. |23|. Some first calculations related to this rigid 
rotor quantization have been done already, with encouraging results. A second issue 
is, of course, the collective coordinate quantization of the (infinitely many) remaining 
symmetries. This point certainly requires further study. A pragmatic approach could 
assume that a more realistic application to nuclei requires, in any case, the inclusion 
of more interactions (even if they are in some sense small), breaking thereby the huge 
symmetry explicitly. Nevertheless, the quantization of the volume-preserving diffeo- 
morphisms may be of some independent interest, although the solution of this problem 
might be difficult. Finally, the semi-classical quantization of the remaining degrees of 
freedom, which are not symmetries, probably just amounts to a renormalization of the 
coupling constants in the effective field theory. These are usually taken into account 
implicitly by fitting the coupling constants to experimentally measured quantities. 
In any case, we think that we have identified and solved an important submodel in the 
space of Skyrme-type effective field theories, which is singled out both by its capac- 
ity to reproduce qualitative properties of the liquid droplet approximation of nuclei, at 
least at the classical level, and by its unique mathematical structure. The model directly 
relates the nuclear mass to the topological charge, and it naturally provides both a finite 
size for the nuclei and the liquid droplet behaviour, which probably is not easy to get 
from an effective field theory. So our model solves a conceptual problem by explic- 
itly deriving said properties from a (simple and solvable) effective field theory. Last 
not least, our exact solutions might provide a calibration for the demanding numerical 
computations in physical applications of more generalized Skyrme models. 
Given this success, it is appropriate to discuss the circumstances which make the model 
relevant. First of all, from a fundamental QCD point of view, there is no reason to ne- 
glect the sextic term, just like there is no reason to ignore the quadratic and quartic 
terms L2 and L4. So the good properties of the L06 model seem to indicate that in cer- 
tain circumstances the sextic term could be more important than the terms L2 and L4. 
The quadratic term is kinetic in nature, whereas the quartic term provides, as a lead- 
ing behaviour, two-body interactions. On the other hand, the sextic term is essentially 
topological in nature, being the square of the topological current (baryon current). So 
in circumstances where our model is successful this seems to indicate that a collective 
(topological) contribution to the nucleus is more important than kinetic or two-body 
interaction contributions. This behaviour is, in fact, not so surprising for a system at 
strong coupling (or for a strongly non-linear system). A detailed study of the general- 
izations mentioned above, or of the more conceptual cosiderations of this paragraph, is 
beyond the scope of this letter and will be presented in future publications. 
Finally, let us briefly mention a recent paper [24|, which appeared after completion of 
this letter. There, a generalized Skyrme model saturating a Bogomolny bound is de- 
rived along completely different lines. The model of that paper consists of a Skyrme 
field coupled to an infinite tower of vector mesons, where these vector mesons may 
be interpreted as the expansion coefficients in a basis of eigenfunctions along a fourth 
spatial direction. Simple Yang-Mills theory in four Euclidean dimensions is the mas- 
ter theory which gives rise to the generalized Skyrme model via the expansion into the 
eigenfunctions along the fourth direction, and the Bogomolny equation for the latter is 
a simple consequence of the self-duality equations for instantons in the former theory. 
If only a finite number of vector mesons is kept, the topological bound is no longer 
saturated, but already for just the first vector meson, the energies are quite close to 
their topological bounds. This latter observation might be in some sense related to the 
results for our model, because integrating out the vector meson produces precisely the 
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sextic Skyrme term in lowest order. One wonders whether it is possible to integrate out 
all the vector mesons, which should lead directly to a topological (Bogomolny) version 
of the Skyrme model. 
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